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(57) Abstract: Material properties such as stress in a ferromagnetic material may be measured using an electromagnetic probe. 
While generating al alternating magnetic field in the object, and sensing the resulting magnetic field with a sensor, the signals from 
the magnetic sensor may be resolved into in-phase and guadrature components. The signals are affected by both geomaterical pa- 
rameters such as lift-off any by material properties, but these influences may be separated by mapping the in-phase and quadrature 
components directly into material property and lift-off components, and hence a material property and/or the lift-off may be deter- 
mined. The mapping may be represented in the impedance plane as two sets of contours representing signal variation with lift-off 
(A) (for different values of stress) and signal variation with stress (B) (for different values of lift-off), the contours of both sets (A, 
B) being curved. The stress contours (B) at a constant angle. Hence calibration measurements taken along a few contours of each 
set enable 
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Measurement with a Magnetic Field 

This invention relates to a method and apparatus in 
which an electromagnetic probe is used to measure 
5 material properties in a ferromagnetic material, for 

example stress, or to measure the separation of the probe 
from the surface of such a material . 

The stresses in structures such as rails, bridges 

10 and pipelines, complex mechanisms such as vehicles and 
machinery, or simple devices such as struts, cables or 
bearings arise from various causes including changes of 
temperature, and the loads and pressures due to use. 
There may also be residual stresses arising from the 

15 fabrication of the structure or device, and any bending 
that the structure or device was subjected to during 
construction; the residual stresses arising from 
fabrication will also be affected by any stress-relieving 
heat treatment. In some situations (such as pipelines) 

2 0 the principal stress directions can be expected to be in 
particular directions (circumferential and longitudinal) , 
whereas in other situations the stress directions are 
also unknown. A variety of magnetic techniques are known 
to have some sensitivity to stress, although magnetic 

25 measurements are usually also affected by other material 
properties such as microstructure . A way of measuring 
stress in a steel plate is described in GB 2 278 450, 
this method using a probe containing an electromagnetic 
core to generate an alternating magnetic field in the 

30 plate, and then combining measurements from two sensors, 
one being a measure of stress -induced magnetic 
anisotropy, and the other being a measure of directional 
effective permeability (DEP) . The probe is gradually 
turned around so the magnetic field has a plurality of 

35 different orientations in the plate, and these 

measurements are taken at each such orientation. The DEP 
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signals are affected not only by stress, but also by the 
lift-off from the surface (i.e. the gap or separation 
between the probe and the surface) , and so must be 
corrected for lift-off. 

5 

According to the present invention there is provided 
an apparatus comprises at least one probe for measuring 
material properties in an object of ferromagnetic 
material and/or measuring the lift-off of the probe from 

10 the surface of the object, the or each probe comprising 
an electromagnet means, means to generate an alternating 
magnetic field in the electromagnet means and 
consequently in the object, and a magnetic sensor 
arranged to sense a magnetic field due to the 

15 electromagnet means; analysis means for resolving signals 
from the magnetic sensor into a first component and a 
second component which are orthogonal in phase; means for 
mapping the first and second components directly into 
material property and lift-off components; and means for 

20 deducing a material property from the material property 
component so determined, and/or a lift-off value from the 
lift-off component so determined. 

The mapping requires a preliminary calibration, with 
25 a specimen of the material, to determine how the first 
and second components of the signal vary with lift-off 
(at a constant stress) and vary with stress (at a 
constant lift-off) , and deducing from the calibration 
measurements the applicable mapping for any stress and 
30 any lift-off. The signals from the sensor are at the 
frequency of the alternating field, and the components 
may be the components in phase with the current supplied 
to the electromagnet means, and the component in 
quadrature to that. The mapping may be represented in 
35 the impedance plane (i.e. on a graph of quadrature 
component against in-phase component) as two sets of 



WO 03/034054 PCT/GB02/04531 

- 3 - 

contours representing signal variation with lift-off (for 
different values of stress) and signal variation with 
stress (for different values of lift-off) , the contours 
of both sets being curved. The contours or lines of one 
5 set intersect those of the other set at non- orthogonal 
angles. The angles at which the contours for constant 
lift-off (varying stress) intersect any one contour for 
constant stress (varying lift-off) are constant along 
that contour. However the angles of intersection of 
10 different lift-off lines with any one stress line are 

different: the angle of intersection at a fixed lift-off 
varies slightly with stress. Hence measurements taken 
along a few contours of each set enable the positions of 
the other contours of each set to be determined. 

15 

Surprisingly this simple mapping has been found to 
give an accurate representation of the variation of the 
signals with material property (e.g. stress or 
micros tructure) ; more surprisingly it enables these 

20 variations to.be distinguished unambiguously from 

variations arising from lift-off or other geometrical 
variations such as surface texture or curvature. It may 
perhaps be presumed that material property changes cause 
changes in the permeability of the steel in the magnetic 

25 circuit and these have both inductive and resistive 

(lossy) components, whereas geometrical changes such as 
lift-off change the amount of air in the magnetic 
circuit, in which energy dissipation by eddy currents 
cannot occur, so this of itself would by purely inductive 

30 (non-lossy) . Where the object has a coating of a non- 
ferromagnetic material (such as a paint, or a ceramic), 
this coating will separate the probe from the surface, so 
the lift-off measurement may indicate its thickness. 

35 Preferably the electromagnet means comprises an 

electromagnetic core and two spaced apart electromagnetic 
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poles, and the magnetic sensor is preferably arranged to 
sense the reluctance of that part of the magnetic circuit 
between the poles of the electromagnet means. It is also 
desirable to arrange for such measurements to be taken 
5 with a plurality of different orientations of the 

magnetic field, at a single location on the object. This 
may be achieved using a single probe that is rotated at 
that location, measurements being taken with different 
orientations of the probe, or using an array of probes of 

10 different orientations that are successively moved to 
that location. In either case, the sensor or sensors 
provide a measure of the permeability of the material 
through which the flux passes between the poles, and so 
provide a signal indicative of the effective permeability 

15 of the material; the corresponding measurements at 

different probe orientations at a location on the object 
hence indicate the effective permeability in different 
directions. The signals from this sensor may be referred 
to as a reluctance, or 1 flux linkage 1 , signal. 

20 

The probe, or at least some of the probes, may also 
include a second magnetic sensor between the two poles 
and arranged to sense magnetic flux density perpendicular 
to the direction of the free space magnetic field between 

25 the poles. This second sensor would detect no signal if 
the material, were exactly isotropic; however stress 
induces anisotropy into the magnetic properties of the 
material, and so the signals received by the second 
sensor are a measure of this stress -induced magnetic 

30 anisotropy or 'flux rotation'. The variations in the 

flux rotation signals at different probe orientations, at 
a location on the object, enable the directions of the 
principal stress axes to be accurately determined. The 
flux rotation signals can also be related to the stress. 

35 

The flux linkage signal from the or each probe is 
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preferably backed-off, i.e. processed by first 
subtracting a signal equal to the signal from that sensor 
with the probe adjacent to a stress-free location. The 
backed-off signal is then amplified so the small changes 
5 in the flux linkage signal due to stress are easier to 
detect. This backing off is performed after resolving 
into in-phase and quadrature components but before 
performing the mapping. Preferably the signals from the 
or each probe are digitized initially, and the backing- 
10 off and resolution are performed by analysis of the 
digital signals. 

To achieve penetration below the surface of the 
ferromagnetic object it is desirable to operate at 

15 alternating frequencies less than 200 Hz, for example 
between 5 Hz and 100 Hz (which in mild steel provide 
penetrations of about 5 mm and 1 mm respectively) . In 
other situations, where such depth of penetration is not 
required, higher frequencies can be used, for example up 

2 0 to 150 kHz for a penetration of only about 15 nm. The 
depth of penetration may be represented by the skin 

depth, 5 = 1/V( TtHoHrfk) , where ^ 0 is the permeability of 

free space, \x r is the relative permeability of the 
material, k is its electrical conductivity, and f is the 
25 frequency. The frequency should be such that the skin 
depth is much less than the thickness of the object. 

Generally, the more different probe orientations are 
used for taking measurements the more accurate the 
30 determination of stress levels and principal axes can be. 
In many cases the principal stress axes can be assumed 
to be aligned in particular directions - axial and 
circumferential directions in the case of a pipe, for 
example - so that the signal maxima for flux linkage 
35 signals would be expected to be along these directions, 
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and the signal maximum for flux rotation signals would be 
along the bisection angles between these directions. 

The probe, or at least some of the probes, may also 
5 include a third magnetic sensor (a 1 f lux leakage* sensor) 
between the poles and arranged to sense magnetic flux 
density parallel to the free space magnetic field. This 
third sensor will detect any flux leakage, this being 
influenced by changes in material properties, lift-off, 
10 and cracks. As with the flux-linkage sensor, 

measurements at a location are preferably made at 
different probe orientations. 

The invention will now be further and more 
15 particularly described, by way of example only, and with 
reference to the accompanying drawings, in which: 

Figure 1 shows a diagrammatic view of an apparatus 
for measuring stress or lift-off; • 

20 

Figure 2 shows a longitudinal sectional view of a 
probe for use in the apparatus of figure 1; 

Figure 3 shows graphically the variation of the 
25 backed-off quadrature and in-phase components of the flux 
linkage signal with variations of lift-off, and with 
variations of stress. 

Referring to figure 1, an apparatus 10 for measuring 
30 stress and/or lift-off includes a sensor probe 12 
comprising sensors for both flux linkage and flux 
rotation, the probe 12 being attached to an electric 
motor 14 which can be held by an operator, so the motor 
14 can turn the probe 12 with one end adjacent to a 
35 surface of a steel object 16 in which the stress is to be 
determined. The sensor probe 12 and motor 14 are 
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connected by a 2 m long umbilical cable 17 to a signal 
conditioning/probe driver unit 18. The unit 18 is 
connected by a long umbilical cable 19 (which may for 
example be up to 300 m long) to an interface unit within 
5 a microcomputer 20, which has a keyboard 21. Operation 
of the apparatus 10 is controlled by software in the 
microcomputer 20 . 

The interface unit within the microcomputer 2 0 

10 generates sine and cosine functions at an angular 

frequency selectable by software, and buffers the sine 
waveform for transmission to the unit 18 for driving the 
probe 12. The amplitude of the transmitted waveform is 
also selectable by software. It also provides signals to 

15 control the motor 14 and hence the angular position of 
the probe 12. The interface unit also provides control 
signals to the unit 18 to select which of the signals 
available from the probe 12 is to the transmitted for 
analysis. It demodulates the selected input signal (flux 

20 linkage or flux rotation) to derive its in-phase and 

quadrature components, filters the demodulated signal to 
remove high frequency components and to reduce noise, and 
converts the analogue signals to digital form for input 
to the computer 20. It also detects the angular position 

25 of the probe 12 from signals provided by a position 
. encoder (not shown) on the motor 14 . 

The long umbilical cable 19 incorporates a coaxial 
cable to transmit the selected signal, and wires to 

3 0 control which signal is selected, to control the motor 
14, to transmit signals from the position encoder, to 
transmit the sinusoidal waveform, and to convey 
electrical power. The unit 18 converts the drive 
waveform from a voltage to a current drive for the probe 

35 12; buffers and amplifies the signals from the probe 12; 
and selects which signal is to be transmitted to the 
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microcomputer 20. It also buffers the signals from the 

position encoder for transmission, and drives the motor 
14 in response to control signals. 

5 Referring now to figure 2, the probe 12 is shown 

detached from the motor 14, in longitudinal section 
although with the internal components shown in elevation 
(the connecting wires within the probe 12 are not shown) . 
The probe 12 comprises a cylindrical brass casing 24 of 

10 external diameter 16.5 mm and of overall height 60 mm, 
the upper half being of reduced diameter whereby the 
probe 12 is attached to the motor 14. The upper half of 
the casing 24 encloses a head amplifier 25. The lower 
half encloses a U-core 26 of laminated mu-metal (a high 

15 permeability nickel/iron/copper alloy) whose poles 28 are 
separated by a gap 7.5 mm wide, and are each of width 2,5 
mm, and of thickness 10 mm. (out of. the plane. of the 
figure) . The poles 28 are in the plane of the lower end 
of the casing 24, which is open. Around the upper end of 

20 the U-core 26 is a former on which are wound two 

superimposed coils 30. One coil 30a (which has 200 
turns) is supplied with the sinusoidal drive current from 
the unit 18; the other coil 30b (which has 70 turns) 
provides flux linkage signals. Between the two poles 28 

25 is a former on which is wound a 1670-turn rectangular 
coil 32, about 4 mm high and 6 mm wide, and 6 mm-square 
as seen from below, the windings lying parallel to the 
plane of the figure so the longitudinal axis of the coil 
32 is perpendicular to the line between the centres of 

30 the poles 28. The coil 32 is supported by a support 

plate 34 fixed between the arms of the U-core 26 so the 
lower face of the coil 32 is in the plane of the poles 
28. The coil 32 provides the flux rotation signals. The 
signals from the coils 30b and 32 are amplified by the 

35 head amplifier 25 before transmission to the unit 18. 
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On the same former as the coil 32 there may also be 
wound a 200 -turn rectangular coil whose windings are 
perpendicular to the plane of the figure so the axis of 
the. coil is parallel to the line between the centres of 
5 the poles. This coil would provide flux leakage signals. 
Alternatively such a flux leakage coil may be provided 
somewhat further from the surface, for example above the 
support plate 34. 

10 In operation of the system 10, the motor 14 is 

supported so the lower end of the probe 12 is adjacent to 
the surface of the object 16 and the longitudinal axis of 
the probe 12 is normal to the surface. An alternating 
current of the desired frequency and amplitude is 

15 supplied to the drive coil 30a, so the magnetic field in 
the object 16 oscillates about zero with an amplitude 
much less than saturation. The probe 12 is first placed 
adjacent to a region' of the object 12 where the stresses 
are negligible. The in-phase and quadrature components of 

20 the flux linkage signal (i.e. the component in phase with 
the drive current, and the component at 90° to the drive 
current) received by the microcomputer 20 are each backed 
off to zero, and the backing off values are then fixed. 
During all subsequent measurements the flux linkage 

25 signal components are backed off by these same amounts 
(i.e. subtracting a signal equal to the component 
observed when in a stress-free location) . 

Measurements can be taken by placing the probe 12 
30 adjacent to a region in which material properties such as 
stress are to be measured. The orientation of the line 
joining the centres of the poles 28 (referred to as the 
orientation of the probe 12) is noted relative to a fixed 
direction on the surface. The motor 14 is then energized 
35 to rotate the probe 12, for example in a step-wise 

fashion 10° at a time through a total angle of 360°. At 
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each orientation of the probe 12 the quadrature flux 
rotation signal is measured, and the flux linkage 
components are measured (and backed off) . These 
measurements may be made at different frequencies, for 
5 example with a drive current frequency of 68 Hz when 
measuring flux rotation, and at a frequency of 300 Hz 
when measuring flux linkage. More typically the drive 
current frequency would be the same when measuring both 
parameters . 

10 

It will be appreciated that the signal analysis 
procedure of the invention is applicable with many 
different probes . The probe 12 might for example be 
modified by using a U-core 26 of a different material 
15 such as silicon iron (which can provide higher magnetic 
fields) , or indeed the drive coil might be air-cored. 

The probe might be of a different shape or size, for 
example for inspecting surf ace . stress in a small bearing 
2 0 it may be appropriate to use a probe of diameter as small 
as 3 mm, and to operate at a higher frequency such as 100 
kHz, while for inspecting internal stresses in a large 
steel pipe it may be appropriate to use a probe of 
diameter say 75 mm. 

25 

The flux rotation signals vary sinusoidally with 
probe orientation, so the orientation at which they have 
their maxima and minima can be determined. The 
directions midway between these two orientations are the 

30 directions of the principal stress axes. Measurements of 
the flux rotation signals are therefore useful if the 
principal stress directions are unknown. The flux linkage 
and flux leakage signals also vary sinusoidally with 
probe orientation (in antiphase with each other) and the 

35 values are observed at the principal stress directions. 
If the principal stress directions are already known, 
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then the probe 12 might instead be merely oriented to 
those directions, and the flux linkage measurements made; 
no rotation of the probe 12 would be necessary. 

5 The values of the stresses in the directions of the 

principal stress axes can be determined from the 
experimental measurements of the flux linkage signals 
with the probe' 12 oriented in those directions. This 
requires calibration of the apparatus 10, taking 

10 measurements on a sample of material of the same type as 
that of the object 16, while subjecting it to a variety 
of different stresses. This may be done with a 
rectangular strip sample in a test rig, measurements 
being made at the centre of the sample where the 

15 principal stress direction is aligned with the axis of 
the test rig. Referring to figure 3 this shows the 
backed-of f flux linkage in-phase and quadrature 
components obtained in such a test rig, the measurements 
being made with a drive frequency of 70 Hz, and the > 

20 specimen being a steel bar. A first set of measurements 
were made at progressively larger values of lift-off, L, 
but with no stress, S. This gives the lift-off line or 
contour A, the lift-off varying between 0 and 220 jim. 
Similar lift-off lines A are obtained for other fixed 

25 values of stress, those for S = 250 MPa tension and 

compression being shown. Measurements were then made at 
a range of different fixed values of lift-off, L, with 
varying stresses, S (both compression and tension), 
providing the contours B. 

30 

It will be appreciated that the contours A are 
curved, and the contours B are not orthogonal to the 
contours A, but that they intersect at angles that are 
constant along any one lift-off line A. Consequently it 
35 is only necessary to make calibration measurements 

sufficient to plot one such contour B and two or three 
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such contours A, and the shapes of the other contours can 
be predicted. It has also been found that the 
intersection angles, and the curvature of the contours, 
are characteristic of the material. 

5 

After calibrating the probe 12 in this manner, 
measurements of stress or of lift-off can be readily made 
from observations of flux linkage signals (resolved and 
backed off) , as the contours enable the changes due to 

10 lift-off to be readily distinguished from changes due to 
stress. Any particular position in the impedance plane 
(i.e. in the graph of quadrature against in-phase 
components) corresponds to a particular value of stress 
and a particular value of lift-off. The mapping between 

15 (in-phase, quadrature) coordinates and (stress, lift-off) 
coordinates, may be carried out graphically, referring to 
such contours, or by calculation. For example if the 
flux linkage .signal has the in-phase and quadrature 
components of the position marked X,;. this corresponds to 

20 a lift-off of about 80 jam and a stress of about -125 MPa. 
Alternatively the value X may be translated along a 
contour B of constant lift-off to find the two components 
(at position R) for zero stress, and so to find the lift- 
off. Similarly this value X may be translated (along the 

25 broken line Y) along a contour A of constant stress to 

find the in-phase and quadrature components at position Z 
for zero lift-off. 

Determination of the positions R and 2 may be 
30 carried out graphically, as shown in the figure. 

Another method is to deduce theoretical equations for how 
the signal will vary with lift-off and material 
properties, taking into account the reluctance of the 
three parts of the magnetic circuit: the U-core 26, the 
35 air gaps, and the magnetic path through the object 16. 
To provide an accurate fit to the experimental data, it 
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would be desirable also to take into account any losses 
in the core 26 and hysteresis losses in the object 16, 
and also flux leakage. Alternatively, simple 
polynomials may be used to model the experimental 
5 measurements; although this approach provides little 

insight into the physical phenomena it has been found to 
be the most accurate approach. This latter approach may 
be modified, based on the experimental observation that 
the intersection angles along any lift-off line and a set 
10 of stress lines remains constant. 

The contours of constant stress, a (referred to as 
lift-off lines) and contours of constant lift-off, L 
(referred to as material or stress lines) constitute a 
15 non-linear two-dimensional matrix on the impedance plane. 
The requirement is to determine for any subsequent 
impedance measurement the independent vector lengths on 
the impedance plane corresponding to a* change in material 
and a change in lift-off. That is to say, the 

20 measurement, (x m ,y m ) corresponding to the position X must 
be separated into the two specific components (x' a ,y' a ) 

at position Z, and (x'L#y'L) at position R, where these 
two components must lie on the zero lift-off stress line 
and the zero stress lift-off line respectively. 

25 

This analysis may be achieved in the following 
stages: 

(A) General analysis for the lift-off /stress matrix 
3 0 measurements 



(1) The lift-off /stress matrix is parameterised using 
simple polynomial equations. Typically this matrix may 
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consist of between 3 and 5 lift-off lines and between 5 
and 10 stress lines. 

(2) All the intercepts between the lift-off and stress 
5 lines are then determined. 

(3) Cubic fits are carried out for both in-phase and 
quadratic intercept values against the actual lift-off 
levels used when measuring the original stress lines 

10 (provided there are at least 4 measured stress lines) . 

(4) Stress lines can then be evaluated for any arbitrary- 
lift-off using the cubic fit coefficients (provided there 
are at least 3 measured lift-off lines) . 

15 

(B) Specific analysis for each subsequent measurement X 

(5) A new lift-off line is evaluated which passes 
exactly through the measurement point. 

20 

(6) The intersection, Z, of this new lift-off line with 
the zero lift-off stress line is determined, (x' a ,y' a ) . 

(7) Find the stress line that passes exactly through the 
25 measurement point. 

(8) The intersection, R, of this new stress line with 
the zero stress lift-off line is determined, (x'l/Y'l) • 

30 (9) Alternatively if stages (3) and (4) have been 

carried out then the new stress line can be compared to 
the cubic fit coefficients to determine the absolute 
lift-off level. 
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The above approach will now be described in greater 
detail, using the numbering of the steps given above. 

(1) All the lift-off (xl/Yl) and stress (x a ,y a ) lines 
5 are first fitted with quadratic equations: 



y i (*) = a L [a)x 2 L + b L {a)x L + c L (a) 



(l) 



10 (where the coefficients are different for each different 
value of a) , and 



y.{L)= a^Ly a ^b 0 {L)x^c a {L) 



(2) 



15 (where the coefficients are different for each different 
value of lift-off) 

It' can be useful to evaluate the intersection angle 
between the lift-off and stress lines, y, at an arbitrary 
20 intersection position (x',y'): 



tan / = 



1 + {2a a {L)x>+b a (L)\2a L (o)x'+b L (a)) 



(3) 



For the lift-off lines we shall specifically 
25 identify the most extreme cases that mark the edges of 

the matrix, for example those for compressive and tensile 
stresses of 250 MPa in the figure. 



30 



where a c is always the largest compression, and 



(4) 
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Vl ) = *l fa )xl + b L (a t )x L + c L (a, ) ( 5 ) 

where at is the largest tensile stress (or at any rate 

the tensile stress giving the largest difference in 
5 signal from that at zero stress) . Also the lift-off line 
measured at zero applied stress, 

yM=a L (0)xl+b L (0)x L+ c L {0) (6) 

10 For the stress lines we shall identify the one at zero 
lift-off, 

y*(0)=a a {0y a +b a {0)x„+cM 

15 (2) The intercepts between measured lift-off and stress 
lines are. determined by setting ( xl , y l ) and (x a ,y CT ) equal 

to (x CT L,y a L) in equations (1) and (2), and solving using 
the fitted polynomial coefficients: 



(b>)-b L (a))±>^o)-b^^ 
20 X ° L= 2(a o( L)-a L (o)) 



then y G h from (1) or (2) . (9) 

Care must be taken to select the solution of equation (8) 
25 which lies within the bounds of the lift-off /stress 
matrix. 

(3) The intercepts, (x oL , Yah) f along each lift-off line 

are fitted to cubic equations to parameterise the true 
30 lift-off values using the 5 to 10 experimental lift-off 
values . 
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to. = P*W + q*W + r X°Y + s x {<r) 



(10) 



y+=P,W +Vy(<7)L 2+r y(°) L + Sy(°) 



(11) 



5 

where the two sets of coefficients p, q and r for x and y 
are numbers that are different for each value of stress 
(three in this example) . 

10 (4) A general complete set of stress lines can then be 
determined using the coefficients in equations (10) and 
(11) for all lift-offs within the measurement range and 
for small extrapolations beyond that range. I.e. 

evaluate general (x a L,y a L) using equations (10) and (11) 
15 for each measured lift-off line, then for each lift-off, 

L=L' , collect all intercepts (x a i,i , y aL . ) and fit 

quadratics to these as in equation (2) but now we have 
generated stress lines for all specified lift-off values, 
these lines making smooth transitions between the 
20 measured stress lines. 

(B) Specific analysis for each measurement on component 
or plant . 

25 (5) For a measurement (x m ,y m ), the lift-off line which 

passes directly through it can be expressed as (referring 
to equations (4) - (6) ) : 



y L = k (0) + (a, ) - a L (a c ))>, 2 + 

[b L (0) + A(b L {a t ) - b L {cr c ))K + k (0) + A {c L (<x, ) - c L (<r c ))] 



(12) 



WO 03/034054 PCT/GB02/04531 

- 18 - 



where the parameter, A, represents the relative position 
across the lift-off /stress matrix along the stress lines, 
and is given by: 



_ y m -a L (0)x 2 m - b L (0>c m - c L (0) 

A " [{a L {<J t )-a L {a c ))xl ^^^tkR + M*,)-c £ (<0)] 



In equation (12) the expressions within square brackets 
are numerical coefficients, say Ai, Bi and Ci. This 

parameter A effectively subdivides the area where there 
10 are contours, giving equally spaced lift-off lines with a 
smooth variation, based on just three measured lift-off 
lines. Note that this analysis does not use the 
information determined in steps (2), (3), or (4) above. 

15 (6) The intersection of the specific lift-off line 

(12) with the zero lift-off stress line, (7), (x' 0 ,y' 0 ), 
is given by: 



, (b„(0)-B,)± V(B, -b o (0)) 2 -4(a o (0)-A,Xc o (0)-C,) 

2(a o (0)-A,) 

20 then y' c from equation (7) by substitution. (15) 



Here again care must be taken to select the correct 
solution. 



25 (7-9) Find the stress line that passes exactly through 
the measurement point. In step (4) above, all general 
stress lines (and their coefficients) were determined. 
Therefore it is only necessary to determine which line 

includes the measurement point, (x m ,y m ) . There are 

30 several ways to do this, one being to determine the 

shortest distance, D, from the point to each stress line: 
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D(L)=^ Xm - X *AL)y + {y m -y%(L)y (16) 



where the (x* CT ,y* a ) are the closest corresponding 
positions on each stress line given by the solution to 
5 the simultaneous equations: 

v *(L)-v x m~ x o*( L ) (17) 

(2a 0 (L)x 0 *(L) + b 0 (L))(a 0 (L>c^ *(L)+b 0 (L)x 0 *( L ) + c 0 (L)-y m ) + (x 0 *(l)-x o )=0 
10 (18) 

Thus for the range of L we find £>(£)Ln and so determine 
the absolute lift-off, L, for the measurement. 

15 An advantage of this method, using stress lines 

determined using the coefficients in equations (10) and 
(11) for all different values of lift-off (as outlined in 
step (4) above), is that it makes use of data from all 
the different experimentally-determined stress lines (of 

20 which there are seven shown in the figure, from L= 0 to 

L=220 jam) . 



Alternatively if absolute lift-off value are not 
used and stages (3) and (4) have not been carried out, 
25 then by working with polynomial intersections in a manner 
analogous to that for the lift-off lines, the 
intersection point R can be obtained. 

An alternative analysis approach takes into account 
30 the experimental observation of a constant intersection 
angle, y, along each lift-off line, in determining the 
equations of the lift-off lines and stress lines. The 
result is more robust to extrapolation beyond the 
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calibrated lift-off range and requires fewer measurements 
of stress and lift-off lines to construct. 

The measured lift-off lines are fitted with simple 
5 quadratic equations as described by equations (4) to (7) . 
These are typically taken so as to represent the two 
extreme values of stress and the zero stress lift-off 
line respectively (three lines, although more can be 
used) . Only a single stress line needs to be measured, 

10 and this can be at any lift-off, Lo, equation (2) . The 

intersection angles, 7(a), of the lift-off lines along the 
stress line are evaluated from equation (3) . Let the 

tangent angle to any lift-off line be <xl(ct) and that to 



15 



20 



the stress line be p a (Lo) . Then: 

tana L (a)^2a L (a)x L +b L (a) . (19) 

tan p 0 (L 0 ) - tan(a L (a) - Y (a)) = tanftan" 1 (2a L (a>c L + b L (a)) - Y (a)]= 2a G (l 0 >c l + b c (l 0 ) 

(20) 

Expanding the more complicated tangent in equation (20) : 

la ,(„), + b fr)-tmM 2 aSh)x b M m) 
1 + (2a L [a)x L + b L (crjjtan y[ij) 

25 In equation (21) all the polynomial coefficients and 

the gamma angles are already known. However we know that 
y(a) is independent of the lift-off, and so we can 
generalise equation (21) to all lift-off values, L: 

\ + {2a L [a)x L + b L (crjjtan y(a) 
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where now the stress line coefficients and the intercept 
between this line and any lift-off line (RHS of equation 
(22)) must be established. 

5 To do this we simply solve the simultaneous 

equations (22) using coefficients from 2 or 3 lift-off 
lines together with equations (2), (4), (5) and (6) to 
determine the stress line coefficients and intercept 

positions. Note that any x position, xl, can be chosen 

10 on a lift-off line and the corresponding y position and 
stress line coefficients determined for that case. 

Hence the position Z can be determined, eliminating 
the effect of any lift-off, and this can be related to 

15 stress. The value of stress found in this way is, it 
will be appreciated, the uniaxial stress that would 
provide that value of the flux linkage signal. If the 
stresses are actually biaxial, then a further calibration 
must be carried out with a cross -shaped sample in a test 

20 rig, flux linkage signal measurements being made at the 
centre of the sample where the principal stress 
directions are aligned with the axes of the test rig. 
Hence a graph or map may be obtained for a range of 
values of stress on one principal axis (say the x-axis) 

25 and for a range of values of stress in the other 

principal axis (say the y-axis) , with contours each of 
which shows the values of biaxial stress that give a 
particular value of apparent uniaxial stress along the x- 
axis; and a similar map may be obtained with contours 

30 showing values of biaxial stress that give a particular 
value of apparent uniaxial stress along the y-axis. 
Hence from measurements of apparent uniaxial stress along 
the two principal stress axes obtained as described 
earlier, the biaxial stress can be determined. 
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It will again be appreciated that the biaxial stress 
may be determined either graphically or by calculation in 
this way. Apparent values of uniaxial stress (in MPa) 
may be used for this purpose, or alternatively the 
5 numerical value of the flux linkage signal (in mV) , 
. either the in-phase or quadrature value, obtained by 
eliminating the effect of lift-off as described in 
relation to figure 3, may be used. 

10 For some materials, such as hard steel, the biaxial 

stress contours for different values of apparent uniaxial 
stress (or of flux linkage signal) along the x-axis, and 
those along the y-axis, intersect at large angles over 
the entire biaxial stress plane. Finding the 

15 intersection of two such contours, and so the true value 
of the biaxial stress, can be done easily. However for 
mild steel the two sets of contours (particularly in the 
tensile/tensile quadrant) are almost parallel, so it is 
very difficult to locate their intersection. In this 

20 case measurements with a different stress-dependent 

variable are helpful, as they can provide a third set of 
contours; the flux rotation signals are desirably used 
for this purpose. 

25 It will be appreciated that the invention enables 

lift-off to be distinguished from changes in material 
properties (such as stress) , over a wide range of values 
of both lift-off and material properties, as it does not 
require an assumption that the changes are linear in the 

30 impedance plane. The probe described above is just one 
type of electromagnetic probe for which the invention is 
applicable. The electromagnet means may include a 
ferromagnetic core on which a drive coil is wound, such 
as a C-core or a rod, or may be air cored. The sensor may 

35 be a sensing coil wound on the same core, or a sensing 

coil positioned near to a pole, or may indeed be the same 
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coil as the drive coil. There may also be several sensor 
coils, for example in an array. 

Measurements of lift-off require at least a region 
5 of surface on which there is no coating, either on the 
object 16 or on a test piece of the same material, so 
that a measurement at zero lift-off can be made. 
Subsequently the thickness of any coating of non- 
ferromagnetic material can be measured, as it is equal to 

10 the lift-off if the probe is in contact with the coating. 
As a rule, a probe smaller than that described above 
would be more sensitive to lift-off and less sensitive to 
material properties, so that if the primary intention is 
to measure lift-off then a probe of external diameter say 

15 8 mm or 4 mm might be used. The shape of a surface, for 
example any undulations in the upper surface of a rail 
head, may therefore be measured in a. non-contact fashion 
by scanning a single probe over the surface and 
monitoring the changes of lift-off determined, as 

20 described above; alternatively an array of probes might 
be used for this purpose, and might also be scanned in 
this fashion. When taking measurements on a surface that 
is curved, it is desirable for the face of the probe to 
have a similar curvature. 

25 

It will thus be appreciated that the present 
invention is particularly suited to the determination of 
stress in a ferromagnetic material, as changes in stress 

have a large effect on the relative permeability, \i r , 
30 while having a negligible effect on conductivity. A 
single specimen can be used for generating the 
signal/stress data and the stress/lift-off data, by 
subjecting the specimen to a range of different uniaxial 
or biaxial stresses. The resulting data can subsequently 
35 be used in taking measurements on many different objects 
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of the. same type of material. In general there is a need 
to determine three independent parameters - the direction 
of the principal axes, and the stresses along each 
principal axis. It is hence desirable, in general, to 
5 take measurements with different orientations of the 
probe . 
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Claims 

1. An apparatus comprising at least one probe (12) for 
measuring material properties in an object (16) of 

5 ferromagnetic material and/or measuring the lift-off of 
the probe (12) from the surface of the object (16) , the 
or each probe (12) comprising an electromagnet means 
(26) , means (30a) to generate an alternating magnetic 
field in the electromagnet means (26) and consequently in 

10 the object (16) , and a magnetic sensor (30b) arranged to 
sense a magnetic field due to the electromagnet means; 
analysis means (20) for resolving signals from the 
magnetic sensor into a first component and a second 
component which are orthogonal in phase; and 

15 characterised by means (20) for mapping the first and 
second components directly into material property and 
lift-off components; and means (20) for deducing a 

« material property from the material property component so 
determined, and/or a lift-off value from the lift-off 

2 0 component so determined. 

2. A method using at least one probe (12) for measuring 
material properties in an object (16) of ferromagnetic 
material and/or measuring the lift-off of the probe (12) 

25 from the surface of the object (16) , the or each probe 
(12) comprising an electromagnet means (26) and a 
magnetic sensor (30b) arranged to sense a magnetic field 
due to the electromagnet means, the method comprising 
arranging the probe (12) adjacent to the object (16) and 

30 activating the electromagnet means (26, 30a) so as to 
generate an alternating magnetic field in the 
electromagnet means (26) and consequently in the object 
(16) , resolving signals from the magnetic sensor (30b) 
into a first component and a second component which are 

35 orthogonal in phase; and characterised by mapping the 
first and second components directly into material 



« 
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property and lift-off components; and deducing a material 
property from the material property component so 
determined, and/or a lift-off value from the lift-off 
component so determined. 

5 

3. A method as claimed in claim 2 also comprising 
calibrating, with a specimen of the material, how the 
first and second components of the signal vary with lift- 
off (at a constant stress) and vary with stress (at a 

10 constant lift-off) , and deducing from the calibration 
measurements the applicable mapping for any stress and 
any lift-off. 

4 . A method as claimed in claim 3 wherein the contours 
15 in the impedance plane representing signal variation with 

lift-off (for different values of stress) and the 
contours representing signal variation with stress (for 
different values of lift-off), are curved, the contours 
of. one set intersecting the contours of the other set at 
20 non- orthogonal angles. 

5 . A method as claimed in any one of claims 2 to 4 
wherein the two sets of contours intersect at a constant 
angle along any one contour representing signal variation 

25 with lift-off. 

6 . A method as claimed in any one of claims 2 to 5 
wherein the electromagnet means comprises an 
electromagnetic core (26) and two spaced apart 

30 electromagnetic poles (28), and the magnetic sensor (30b) 
is arranged to sense the reluctance of that part of the 
magnetic circuit between the poles (28) of the 
electromagnet means. 



35 



7. A method as claimed in claim 6 comprising taking 
such measurements on the object with a plurality of 
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different orientations of the magnetic field, at a single 
location on the object. 

8. A method as claimed in claim 4 or claim 5 wherein 

5 the contours are represented by polynomials in performing 
the mapping. 

9. A method as claimed in claim 5 wherein the contours 
are presented by polynomials in performing the mapping, 

10 the constant intersection angle along a contour 

representing signal variation with lift-off being taken 
into account in generating the polynomials. 
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